Critical Current and Topology of the Supercooled Vortex State in NbSe2 
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We study the behavior of the critical current, Ic{H,T), of pure and Fe doped NbSe2 crystals in 
the denominated disordered vortex region, limited by the critical field Hc2{T) and the field Hp{T) 
at which the peak effect in Ic{H,T) is detected. The critical current follows an individual pin- 
ning response as demonstrated by its field independent universal function of the superfiuid density. 
Transport measurements combined with Bitter decorations show no evidence of the existence of an 
amorphous phase in the high temperature region. 
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I. INTRODUCTION 

Although the peak effect in the critical current, 
Ic{H,T), at a field Hp{T) close to the upper critical 
field, Hc2{T), has been observed years ago its origin is 
still subject of debate. The traditional interpretation 
of this effect is based on the explanation suggested by 
Pipparda and on the formal treatment made by Larkin 
and Ovchinnikov (L0)o. In this scenario the peak effect 
appears as a response of the vortex lattice to the pre- 
sence of quenched disorder. The peak occurs due to the 
softening of the elastic constants in a crossover from local 
to non-local elasticity when approaching Hc2{T). Using 
this formalism the behavior of Ic{H,T) at the peak in 
NbSe2 has been explained.El 

The elastic description, of the peak effect has been cri- 
ticized by other authorsu emphasizing the necesary dy- 
namic aspects of the phenomenon. An extensive anaJysis 
through transport measurements in NbSe2 showf do the 
plastic fingerprint predicted in 2D simulationsQ'El in the 
H — T region where the peak effect appears. In Ref. ||a 
phase diagram with a liquid vortex phase for H > Hp{T) 
is proposed. 

Recent electrical transport measurements using 
Corbino-contact configuration have suggestedQ other in- 
terpretation of the peak effect phenomenon. An abrupt 
decrease of IdH, T) at Hp{T) is considered to be the ma- 
nifestation of a genuine first order thermodynamic phase 
transition from a disordered {H > Hp{T)) to an ordered 
vortex state {H < Hp{T)). In this scenario the width 
of the transition detected using the conventional four 
contact configuration is duecl to the injection of a dis- 
ordered sZjOEtex phase through the surface. It is widely 
acceptedoEEl that the critical current is history depen- 
dent for H < Hp{T) while for Hp{T)<H <Hc2{T) it is 
reversible. A dramatic example of historyj-dependence of 
I — V curves has been recently reportedtj using a fast 
transport measurement technique. This technique al- 
lows to define a large critical current when the sample 
is cooled, in the presence of a field (FC), down to tem- 
peratures well below that of the peak in Ic{H,T). On 
the other hand, when using a slow transport technique a 



low critical current is detected. When the field is applied 
after zero field coohng the sample (ZFC) to temperatures 
below that of the peak in the critical current, the same 
low value of Ic{H,T) is obtained independently of the 
technique, fast or slow, used to measure it. 

Within the frame of this work the most relevant pro- 
posed properties of the mixed state associated with the 
presence of the peak effect can be summarized as follows: 
(a) The vortex phase nucleates at Hc2{T) as a topolo- 
gically disordered structure with a large critical current 
that remains in thermodynamic equilibrium, down to the 
temperature where the peak in Ic{H,T) takes place; (b) 
when the sample is FC to lower temperatures the disorde- 
red vortex structure is supercooled as a metastable state 
with a large Ic{H, T); (c) the thermodynamic equilibrium 
state obtained by ZFC experiments at low temperatures 
corresponds to a topologically ordered phase with low 
critical current. 

Independently of the different interpretations of the 
peak effect there is consensus on the existence of an equi- 
librium disordered vortex state at temperatures above 
that of the peak, where transport measurements show 
no history dependence. It has been claimecO'Ll that the 
disordered equilibrium phase shows a reentrance at fields 
low enough. In particular, in Fe doped samples the reen- 
trant region has been found a a field of 1000 Oe.Q 

The history independence of the transport properties 
when H > Hp{T) suggests that the vortex structure is 
displaced homogenously through the sample when the 
measuring current exceeds the critical one. However, 
the analysis of the field and temperature dependence of 
Ic{H,T) in the so called disordered phase is still lack- 
ing. On the other hand, the claim that in FC experi- 
ments the disordered phase remains in a metastable state 
at low temperatures provides us the possibility of using 
the Bitter decoration to visualize its topology. This is 
particularly important because despite the claim of the 
presence of a disordered phase based on the analysis 
Ic{H,T) measurements, there is no trivial relationtiT 
between critical current and topological order. The criti- 
cal current is determined by vortex displacements in the 
small length scale associated with the range of forces of 
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the pinning potential (superconducting coherence length 
^(T)), while the vortex topology responds to a larger 
length scale (lattice parameter), see Ref. H. 

In this work we focus the attention on the study of the 
field and temperature dependence of the critical current 
and its relation with the corresponding vortex structure 
of the proposed disordered phase. This is accomplished 
by complementing the transport measurements with Bit- 
ter decoration of the pristine FC vortex structure before 
any external force is applied. Thus, the Bitter image de- 
picts the vortex structure tested by fast transport mea- 
surements. 



this behavior is difficult to check due to the narrow avail- 
able experimental field range between Hp{T) and Hc-AT). 

In order to verify the behavior predicted by Eq. we 
investigated both, the field as well as the temperature de- 
pendence of Jc{H, T). It should be noticed that Eq. || has 
an explicit dependence in magnetic field when Jc{H,T) 
is plotted as a function of Q. This dependence should 
be evident when measuring the critical current as a func- 
tion of temperature for different fields. In this way it can 
be proved whether the vortex structure corresponds to 
an homogeneous vortex phase described by a correlation 
volume limited in its transverse direction, Rc = aQ. 



II. LARKIN-OVCHINNIKOV THEORETICAL 
BACKGROUND 

Since the LO theory has been shown! to be appropria- 
te to provide a quantitative understanding of IdH, T) in 
the peak effect, we take the theoretical prediction of the 
model as the basic ingredient for the analysis of Ic{H, T). 

This model predicts an enhancement of the critical cu- 
rrent when approaching Hc2{T) due to the softening of 
the vortex lattice associated with the dispersive charac- 
ter of the tilt modulus. The increase of Jc is associa- 
ted with a monotonous reduction of the LO correlation 
volume Vc = R^L^, where Rc is the elastic correlation 
length perpendicular to H and Lc is the parallel one. At 
Hp{T), Rc takes the minimum value compatible with the 
existence of a finite density of vortices, Rc — ao, where 
ao = 1.075($o/S)^^^rjs the vortex lattice constant. In this 
picture it is assumed^ Rc = uq in the whole region of fields 
Hp{T)<H<Hc2{T) OT equivalently Tp{H) <T <Tc2{H). 

In the LO theory the critical force is given by 
BJcVc — fN^/^, where N is the effective number of pin- 
ning centers in Vc and the pinning force, /, is usually as- 
sumed to be proportional to ^P^, where ^ is the Ginzburg- 
Landau order parameter. Close to Hc2{T) the conden- 
sation energy decreases as ^I*^ oc (1 — &)(1 — t) where 
t = T/Tc and b = B/Hc2{T). It is easily seen that 
(l-6)(l-t) = eTc where Q^T~Tc2{H). Within these 
considerations and using the non-local elastic constant 
C|4 = C44(1 -b){l- t\U'^, where CU^ « H"^ is the local 
tilting elastic constantE£l, it is foundElthat for H > Hp{T) 

Vc = Leal = aliAClAlfn"^?'' oc iJ-^^ (i) 

where it has been assumed that close to Hc2{T) all pi- 
nning centers in Vc are effective: N = nVc. Thus, the field 
and temperature dependence of the critical current is 

Jc{H,T) (X H-'^/'^e (2) 

Most of the experiments reported in the literature have 
been made at constant temperature varying field. In this 
case the theory predicts that the critical current decreases 
as (1— &) and the rapid change of Jc{H, T) as a function of 
h has been associated to this field dependence. However, 



III. EXPERIMENTAL 

The experiments were performed in pure and Fe doped 
NbSe2 single crystals. Material parameters for the Fe 
doped sample are Tc = 5.8 K, 5 Tc{lO% - 90%) = 0.2 K, 
^(0) = 78Aand for the pure sample Tc = 7.02 K, 
(5rc(10% - 90%) = O.IK, ^(0) = 89 A. 

The results for the pure and impure samples are quali- 
tatively similar. However, measurements of the Fe doped 
sample have a relative advantage over the pure one due 
to the wider field region Hp{T) <H < Hc2{T). 

The transport properties were measured using the 
Corbino-contact configuration as shown in the insert of 
Fig. 6 (b). A dc current was injected between the cen- 
tral contact and the four contacts indicated in the fig- 
ure. The distance between and V~ is approximately 
0.3 mm and between V'^ and V~ is 0.3 mm, in both cases 
the sample thickness is around 40 /im. 

The sample was thermally anchored to a sapphire 
glued to a copper sample holder that was thermally 
weakly connected to a liquid helium bath at 4.2 K. A 
heater on the sample holder was used to control the tem- 
perature of the sample by means of a FID controller. 




I[mA] 

FIG. 1. Typical I — V curves for the Fe doped sample using 
the Corbino-contact configuration, where tlie critical current 
is marked with an arrow. The measurements were made at 
different magnetic fields and T = 5.52 K greater than Tp{H). 
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Measurements were made either cooling the sample in 
a field, FC experiments, or cooling the sample in zero 
field and then sweeping the field at a given temperature, 
ZFC. 

The critical current was defined as that inducing a vol- 
tage of 10 nV in the Fe doped sample and of 20 nV in the 
pure sample. Moderate changes in voltage criteria to 
define Ic{H,T) do not modify the conclusions reached in 
the paper. 

The Bitter decoration in pure samples was performed 
at 4.2 K in FC experiments following the procedure des- 
cribed in Ref. The magnetic decoration of the Fe 
doped sample was made at 3 K in order to achieve the 
necessary resolution. The ZFC decorations of pure sam- 
ples were made applying the field after cooling the sample 
down to 4.2 K. 



IV. RESULTS AND DISCUSSION 

Figure |l| shows typical I—V curves for different magne 
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FIG. 2. Critical current curves for different magnetic fields 
as a function of temperature for the Fe doped sample: (a), (b), 
(c) and for the pure sample: (d), (e), (f). The position of the 
peak in Ic{H,T) is indicated by an arrow. Figures (c) and (f) 
show the Ic{H,T) curves at the lowest field where the peak 
is detected for each sample. 



tic fields obtained in the Corbino contact configuration 
for the Fe doped sample in the temperature range T > 
Tp{H). The critical current for each magnetic field is 
marked with an arrow. Similar results are obtained for 
the pure sample. For the same value of I—Ic the voltage 
increases with H , as expected. 

Critical current curves as a function of temperature for 
different magnetic fields are shown in Fig. ^jfor Fe doped 
and pure NbSe2 samples. A well defined peak in IdH, T) 
is found for fields higher than 200 Oe and 25 Oe for the 
Fe doped and pure samples respectively (see Fig. ||(c) 
and (f)). Bellow those fields only a change in the slope 
of Ic{H,T) is detected when approaching Tc2{H). We 
have defined Hc2{T) (or Tc2{H)) as the magnetic field 
(or temperature) where the critical current extrapolates 
to zero. The results are shown in Fig. |^ (a) and (b). 

The critical current of both materials obtained with 
the described experimental setup is history independent 
and consequently Ic {H, T) is a single- valued function of 
H and T. An example is shown in Fig. ^ for the Fe 
doped sample. Most of the data reported in this work 
was taken sweeping temperature. 

To verify whether the LO theory (Eq. ||) describes the 
field and temperature dependence of the results shown in 
Fig. 1^, the critical current for different fields is plotted 
as a function of the variable O, see Fig. ^. Contrary 
to what is expected from Eq. |^, the critical current in 
the variable O becomes field independent in a region of 
temperatures higher than that where the peak in lc{H, T) 
takes place : —0.1 < < for the Fe doped sample and 
—0.03 < O < for the pure case. This is clearly shown in 
Fig. |where/c(i?,e) as a function of magnetic field for 
different values of O is plotted. This field independent 
regime for the critical current is obeyed even for the range 
of fields where the peak effect is not observed, see Fig. 
1(a). 
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FIG. 3. Field dependence of Ia{H) a.t T = 5.3 K in the 
Fe doped sample. The filled square symbols depict the peak 
effect for FC measurements while the open circles are for the 
ZFC case. 
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As a consequence of the result described in the pre- 
vious paragraph it is concluded that the LO theory (see 
Eq.^ is not applicable in the O region where the criti- 
cal current is field independent. Besides this, the field 
dependence of Ic{H,Q) at the peak docs not follow the 
behavior predicted by the theory in the limit Rc = a^. 
The predicted increase of Ic{H, 9) at the peak when de- 
creasing the magnetic field is not observed in any of the 
samples. Thus, the behavior of the critical current in the 
peak effect region is not qualitatively described by the 
LO collective theory under the assumption Rc — aQ. 

The results in Figs. ^ and ^ make evident that the 
field and temperature dependence of the critical current 
are only determined by the change of superfluid density 
as described by A,; 2 (r,iJ) oc -^"^ (x (1 - - 6) = 9 Tc, 
when approaching Tc2{H). The lack of explicit field 
dependence of Ic{H,Q) indicates the irrelevance of the 
vortex-vortex interaction to determine the pinning en- 
ergy, suggesting that the critical current is fully deter- 
mined by the pinning of individual vortices: single vortex 
limit. This seems to be a counter- intuitive and surprising 
result in a field region where the distance between cores 
is smaller than the interaction characteristic length. 



Xh{T, H). We will come back to this point in the conclu- 
sions. 

Figure ^ summarizes the results described up to now in 
a B—T phase diagram. In Fig. ^ (a) the line Hp{T) indi- 
cates the field where the peak in Ic{H, T) takes place and 
the line Hs{T) marks the crossover to the single vortex 
regime (lack of field dependence of Ic{H^ 9)) for the Fe 
doped sample. The Fig. || (b) shows the results for the 
pure sample, in this case Hs{T) coincides with Hp{T). 
The new individual vortex regime line Hg (T) is found to 
follow the condition Qs=Ts{H) — Tc2{H) — constant. 

In order to detect a possible reentrance of the line 
Hp(T) in the phase diagram of the Fe doped sample as 
reported in Ref. 0, we have made careful measurements of 
critical current as a function of temperature or magnetic 
field keeping constant field or temperature respectively. 
No evidence of a reentrant Hp{T) was found for the inves- 
tigated sample (the paths followed to detect the possible 
reentrance are indicated with dotted lines in Fig. ^(a)). 
The presence of a peak effect down to 200 Oe and the 
absence of reentrance in the Fe doped sample is a diffe- 
rent behavior than that found for similar samples in Ref. 
(comparable Tc and ^(0)) where the reentrance of the 
high temperature phase takes place at iJ = 1000Oe. 
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FIG. 4. Critical current as a function of the temperature 
(see text) for different magnetic fields, (a) Fe doped sam- 
ple: Dark symbols correspond to the results where no peak in 
Ic{H,Q) is detected. The insert shows the same curves as a 
function of the absolute temperature, (b) Pure sample. The 
data shows no agreement with the LO theory, see text. 



FIG. 5. Plot of Ic{H, 0) as a function of H for different 
0. For the Fe doped sample (a) the single pinning regime is 
found for values of > —0.1 and for the pure sample (b) for 
0>-O.O3. 
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FIG. 6. B — T phase diagrams, (a) Fe doped sample 
{Tc = 5.8K). Hc2{T) is the field where 7c goes to zero. The 
Hp{T) line is the position of the peak of Ic{T,H). The 
line Hs{T) corresponds to a crossover to a single pinning 
regime. Dotted lines indicate the paths followed to study 
the reentrance of the disordered vortex phase, (b) Pure sam- 
ple {Tc = 7.02 K). In this case the Hs{T) line coincides with 
Hp{T). In the insert is shown the Corbino contact configura- 
tion. 

This clearly shows that the reentrant peak effect depends 
strongly on unknown properties of the sample and conse- 
quently on the ubiquity of the supposedjjiisordered phase 
as has been verified by several authors.E£l 

The results of magnetic decorations performed in FC 
experiments at 60 Oe are shown in Fig. ^(a) for the pure 
sample at 4.2 K and in (b) for Fe doped sample at 3K. 
It is important to remark that in this case the FC deco- 
ration for the Fe doped sample is made at a field below 
the lowest field where the peak is detected. In the pure 
case the FC at 60 Oe crosses the Hs{T) line. An iden- 
tical polycrystalline nature of the FC state is observed 
in both type of materials. Following Ref. |l^, this struc- 
ture corresponds to the so called metastable disordered 
phase. The polycrystalline structure is detected indepen- 
dently of the presence of the peak effect at the field where 
decoration is made. This indicates that the FC polycrys- 
talline structure is originated in the region of tempera- 
tures where the single vortex regime dominates. This is 
an interesting result showing that some degree 



FIG. 7. FC vortex structure in NbSe2 observed by mag- 
netic decoration at 60 Oe. Delaunay triangulations where the 
nearest neighbors vortices are bonded and the non-sixfold co- 
ordinated are marked in gray, (a) Pure sample decorated at 
4.2 K and (b) Fe doped sample decorated at 3 K. 

of orientational order is preserved when nucleating the 
vortex structure at Hc2{T) in a single vortex limit. 

The Bitter structure obtained in ZFC experiments at 
4.2 K shows gradients of i? as a result of the presenpj 
of transport currents induced when applying the field.tZI 
Despite this, a polycristalline structure is found in most 
of the sample. The average size of the grains coincides 
with those shown in Fig. M obtained in FC experiments. 



V. CONCLUSIONS 

We have found a universal behavior of the tempera- 
ture and field dependent critical current when approa- 
ching the mean field superconducting transition, Tc2{H). 
The universal character is made evident when plotting 
Ic{H,T) as a function of the temperature difference = 
T - Tc2{H), as shown in Figs. ^ and ^. This shows that 
the superfluid density univocally determines the pinning 
properties of the material in this region of the phase di- 
agram. The lack of explicit field dependence of Ic{H, 8) 
when plotted as a function of O indicates that the vortex- 
vortex interaction plays no significant role in the deter- 
mination of vortex pinning. The results cannot be inter- 
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preted withui the collective pinning theory in the single 
vortex limitB since the distance between vortex cores is 
much smaller than the effective penetration depth, A/j. 
This rather suggests that in this large Xh limit, the in- 
teraction of individual pinning centers with a vortepc line 
overcomes the weak interaction between vortices .113 

It is important to remark that the novel single pinning 
regime is not necessarily associated with the presence of 
the peak effect: despite the peak effect is not detected 
below 200 Oe in the Fe doped sample the crossover to 
single pinning takes place even at lower fields. 

The transport and Bitter experiments make evident 
that the single vortex pinning found at high temperatures 
is compatible with a structure where directional order is 
preserved within the grains. This reinforces the claim 
that the typical length scale associated with critical cur- 
rent is-different from that characterizing the topological 
order .til 

We have not been able to detect a reentrance of Hp (T) , 
neither in the Fe doped sample nor in the pure one. 
Whether this indicates that the phenomenon fades away 
when decreasing field or transforms into a weak signal 
not detected in our samples deserves further research. 

The combined transport and Bitter decoration ex- 
periments together with the results in Ref. |l^ open a 
question mark on the nature of the disordered phase. 
This is of naxticular importance in view of recent data 
interpretedQcj'El in terms of an order-disorder thermo- 
dynamic transition that takes place at the temperature 
where the peak in Ic{H,T) is observed. The results re- 
ported in this paper call for a better definition of what 
is meant by order and disorder characterized by critical 
current measurements. The topology of the FC struc- 
ture as observed by Bitter decoration shows no gradient 
of B indicating the absence of transpjost currents act- 
ing on the bulk of the vortex system.E3E3 Although this 
supports that by FC the high temperature structure is 
supercooled, the results make evident that the high tem- 
perature phase is not amorphous. 

The observed ZFC vortex structure is not an ordered 
equilibrium state; it is consistent with a topological 
unstable-State dominated by the presence of transport 
currents. 113 

In conclusion, through the analysis of the vortex struc- 
ture and the assumption of the existence of the proposed 
supercooling of the high temperature phase in the FC ex- 
periments, no evidence of an order-disorder phase tran- 
sition is found in samples where the transport properties 
show anomalies associated with a first order thermody- 
namic transition!] 
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